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ABSTRACT 

We present mid-infrared observations with the Spitzer Space Telescope of the nearby type II- P 
supernova, SN 2004dj, at epochs of 89 to 129 days. We have obtained the first mid-IR spectra of any 
supernova apart from SN 1987A. A prominent [Nill] A6.64^tm line is observed, from which we deduce 
that the mass of stable nickel must be at least 2.2 x 10~ 4 M . We also observe the red wing of the 
CO-fundamental band. We relate our findings to possible progenitors and favour an evolved star, 
most likely a red supergiant, with a probable initial mass between ~10 and 15 M®. 
Subject headings: supernovae: general — supernovae: individual(SN 2004dj), 



1. INTRODUCTION 

Core-collapse supernovae (SNe) are the end points of 
most stars more massive than ~8Mq. As such, they 
provide a key test of stellar evolution. Furthermore, they 
play a major role in driving the chemical and dynamical 
evolution of galaxies, and have also been proposed to be 
major contributors to dust at epochs when the Uni verse 
was still young (z > 6) (e.g. lTodini fc Ferrarall200l|L 

SN explosions provide unique natural laboratories for 
studying, in real time, the physics of a variety of com- 
bustion, hydrodynamic, nuclear, and atomic processes. 
While SNe constitute important astronomical sources 
over all wavelength ranges, the combination of strong ab- 
sorption by the Earth's atmosphere and the high back- 
ground at mid-IR wavelengths has meant that this re- 
gion has thus far remained inaccessible for the study of 
SNe, apart from the exceptionally nearby SN 1987A (e.g. 
iRoche et al~1ll993t iWooden et al. Ill993|) . 

Nevertheless, the mid-IR region holds the potential of 
providing unique insights into the nature of SN explo- 
sions and the role played by dust in these events. Al- 
though abundance measurements have long been carried 
out using UV/optical spectra, the large number of lines 
coupled with strong Doppler broadening leads to line- 
blending resulting in ambiguities in species identification 
and errors in flux measurement. In contrast, the fewer 
line transitions and much-reduced sensitivity to extinc- 
tion in the mid-IR, allows firm line identifications and 
accurate measurements of intrinsic line strength. Abun- 
dance measurements using fine-structure lines are partic- 
ularly robust as these are largely insensitive to tempera- 
ture. Furthermore, warm dust that may condense in the 
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ejecta emits most strongly in the mid-IR. Moreover, by 
monitoring the mid-IR spectral energy distribution and 
evolution, we may discriminate between pre-existing cir- 
cumstellar dust and newly condensing dust in the ejecta. 

The advent of th e Spitzer Space Telescope (SST; 
I Werner et al. I 120041) with its vastly improved mid- 
IR sensitivity/spatial resolution combination, compared 
with previous instrumentation, has finally opened up the 
possibility of studying typical nearby SNe in the mid-IR. 
In this Letter, we report on the first mid-IR results for 
SN 2004dj, the nearest SN in over a decade. 

SN 2004dj was discovered in the nearby nearly face-on 
spiral galaxy NGC 2403 on 2 004 July 31 by K. It agaki at 
a visual magnitude of +11.2 (Nakano et al. 2 0Q4|). It ap- 
pears within a star cluster ("Sandaee 96: iSandage 1119841 
iMaiz-Apellaniz et alJ I2004D. A spe ctrum obtained on 
2004 August 03 bv lPatat etaTI ll200l revealed SN 2004dj 
to be a type II-P (plateau) SN at an epoch of ~3 weeks 
post-explosion. In w hat follows, we assume a n explosion 
date of 2004 July 10. iKorcakova et alJ 1)200511 present an 
i?-band light curve which shows that the plateau phase 
had ended by about +100 d post-explosion. NGC 2403 
lies within the M81 g roup at a distance of 3.13 Mpc 
ijFreedman et al. l l2001). About a month post-explosion, 
SN 200 4dj was also detected at radio and X-ray wav e- 
lengths (|Stockdale et al.ll2004l IPoolev fc Lewin II2004J) . 

2. OBSERVATIONS 

SN 2004dj was observed with the SST as part of Di- 
rector's Discretionary Time (DDT) Program 226. Here 
we report on the first two of four epochs with all three 
instruments (IRAC, MIPS, IRS) at epochs of +89 to 
+139 d. Two further epochs of IRAC and MIPS images 
were serendipitously obtained from the Spitzer Investiga- 
tion of Nearby Galaxies (SINGS) Legacy project (Spitzer 
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TABLE 1 
Photometry of SN 2004dj 





Epoch t cxp 




Flux (mjy) 
IRAC 








Epoch 


texp 


MIPS 


Date 


(d) (s) 


3.6 


4.5 5.8 


8.0^m 


Date 




(d) 


00 


24/im 


2004 Oct. 07 


+89 150 


9.99±om :( 


7.60+0.03 6.53+0.05 


4.09+0.04 


T2004 Oct. 


12 


+94 


41.9 


1.1+0.3 


t2004 Oct. 10 


+92 240 


10.82+0.03 


8.64+0.03 6.15+0.05 


3.93+0.04 


2004 Oct. 


11 


+96 


165.7 


1.2+0.2 


t2004 Oct. 12 


+94 240 


7.41+0.03 


8.43+0.03 5.72+0.05 


3.39+0.04 


t2004 Oct. 


16 


+98 


41.9 


1.0+0.3 


2004 Nov. 01 


+114 150 


4.56+0.02 


6.67+0.02 4.48+0.04 


2.48+0.04 


2004 Nov. 


06 


+119 


165.7 


1.1+0.2 



Notes. The epochs are for an assumed explosion date of 2004 Jul. 10 IPatat et al.l200*4T) . t eX p is the on-source integration time. Photometry 
was carried out using a 3'.'6 radius aperture for all four IRAC bands. Aperture corrections of 1.12, 1.12, 1.14, and 1.23 were applied to 
channels 1-4 respectively (cf. Table 5.7 of the IRAC data handbook). For the 24/xm MIPS data, we used a 5"6 radius aperture (aperture 
correction = 1.72; MIPS data handbook, Fig. 3.2). Using a 6'.'0 radius aperture for the MIPS 70/mi channel, we estimate a rough upper 
limit of lOmJy. t exp at 70 and 160^tm is, respectively, 125.8 and 25.2 s for our set-up, and 41.9 and 4.2 s for the SINGS set-up. + The errors 
shown above are statistical errors only. Note that systematic errors in the calibration can be as large at 10%, although the relative errors 
are likely to be much smaller, t SINGS data. 
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Fig. 1. — Subsection of a 3.6/^m Spitzer (IRAC) image taken at 
+89 d showing the immediate field around SN 2004dj. Short dashes 
mark the position of the SN (a J2 ooo = 07 ,l 37 m 17+02; <5j 20 oo = 
+65°35'57".8) which lies 10"N and 160"E of the galaxy nucleus. 

PID 159: iKennicutt et al. ll200^ in 2004 October. 

2.1. Photometry 

From the IRAC images it is clear that we have detected 
a source coincident with the position of SN 2004dj (Fig. 

That this source is dominated by emission from the 
SN is evident from its fading with time (Table 0. Aper- 
ture photometry was carried out on the IRAC and MIPS 
(24/xm) post-BCD image mosaics using GAIA. We re- 
peated the procedure for two IRAC epochs (+89 and 
+114 d) using the APEX software applied to the BCD 
frames and found generally consistent results. The re- 
sults of our measurements are given in Table The SN 
was not detected in the 70 and 160 pm. MIPS channels, 
probably due to a combination of high host galaxy back- 
g round and poorer spati al resolution in these channels. 

ISueerman et al.1 l|2005j) have reported the rapid (5-day) 
decline of SN 2004dj at mid-infrared wavelengths. Wc 
point out th at between +89 d and 11 4 d the plateau phase 
was ending ijKorcakova et all 12005(1 which accounts for 
the rapid decline in the 3.6/im and 8.0/im channels dur- 
ing this time. The slower decline at the other two wave- 



lengths is probably due to the emerging CO fundamental 
band emission (see § 13. 1|) . There is also apparent evidence 
of a slower decline at 24 fim, but this only significant at 
the 80% level. We therefore do not discuss it further. 

2.2. Spectroscopy 

Spectra were obtained on 2004 October 24 (+106 d) 
and 2004 November 16 (+129 d) using the Short-Low 
(SL, A = 5.2 - 14.5/iin; R = 64-128) and Long-Low 
(LL, A = 14 - 38^m; R = 64-128) modules of the In- 
frared Spectrograph. The observations were performed 
in staring mode. Total time-on-source for both epochs 
are 609.5s in SL and 629.1s in LL. Module slits were 
positioned relative to a reference star using a moderate- 
accuracy peak-up with the blue (16/im) channel. 

The SL data were preprocessed using version SI 1.0. 2 
of the Spitzer data processing pipeline. Subsequent re- 
duction of these post-BCD data was carried out within 
the FIGARO 4 environment. We removed the back- 
ground emission by differencing the two nod positions. 
For the first order SL spectrum, the sky position (for 
both epochs) unfortunately landed on a cluster of bright 
sources, rendering it unusable. We therefore extracted 
the spectrum from the single uncontaminated nod posi- 
tion. We used the 'tune' tables for the wavelength cal- 
ibration. Flux calibration was carried out in two steps: 
each order was calibrated separately using the FLUX- 
CON keywords provided. We then merged the spectra, 
manually clipping the ends of each order where the noise 
increases significantly. The spectra were scaled to match 
our day +114 IRAC channel 3, 4 photometry. The IRS 
epochs are, respectively, 8d before and 15 d after this 
IRAC epoch, but no other post-plateau IRAC points 
were available. However, the very slow mid-IR evolution 
of SN 1987A around these epochs suggests that the flux- 
ing error introduced by the epoch differences is probably 
small. We applied scaling factors of of x0.84 and x 0.8 to 
the 106 and 129 d spectra respectively. The uncertainty 
in these factors is ~15% which probably dominates the 
overall fluxing error. We repeated the above reduction 
sequence for the +106 d spectrum, but this time starting 
from the BCD data and found excellent agreement with 
the post-BCD results. The spectra are shown in Fig. 
The SN was not detected in the LL modules. 

We supplemented the IRS data with 0.9-2.4 pum spec- 
troscopy carried out on 2004 November 24 (+137d) us- 
ing the LIRIS instrument mounted at the Cassegrain fo- 
cus of the William Herschel Telescope. These data were 
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Fig. 2.— IRS spectra of SN 2004dj. The 106 d spectrum has 
been shifted by +2,5 mjy. The lowermost curve is the 260 d sp ec- 
trum of SN 1987A IMeikle et al.lll98gt I Wooden et~ai~l 11993) : it 
has been arbitrarily scaled to highlight the striking similarity with 
SN 2004dj. The 114d IRAC photometry is also shown. The inset 
zooms in on the 6.2-7.8/im region of the SN 2004dj spectra (the 
129 d spectrum has been shifted by +2.5mjy) showing the profiles 
of the [Nill] line and the emergence of the blend near 7.5/im. 

obtained using the low-resolution ZJ and HK grisms 
in the standard ABBA pattern with 10" nods. The 
data were reduced in the usual fashion; the F6 dwarf, 
BS 3028, was our chosen flux standard. Final fluxing 
was achieved by comparison with field stars in the J- 
and ii-band acquisition images, calibrated using 2MASS 
data. The derived SN magnitudes are J = +f3.62±0.04 
and H = +13. 30 ±0.04. We used these values to derive 
scaling factors {ZJ = xf .f4 and HK = x2.68) which we 
applied to the near-IR spectrum shown in Fig. |31 (The 
2MASS pre-explosion magnitudes of the underlying clus- 
ter are J = +16.04 ± 0.11 and H = +15.74 + 0.12). 

3. ANALYSIS 
3.1. Carbon Monoxide 

The two IRS spectra show little difference in overall 
appearance. The blue- most region (see Fig. [21, is domi- 
nated by a rapidly rising slope which we identify with the 
red wing of the carbon monoxide fundamental (4.65 (an). 
This identification is reinforced by the clear detection of 
(i) the first overtone of CO at ~ 2.3^m (Fig. and (ii) 
the growth with time of the excess flux in channels 2 and 
3 compared with channels 1 and 4 (see Table 0. 

Molecule formation provides a sensitive diagnostic of 
the conditions and degree of mixing in the SN ejecta. 
Cooling by CO sets the temperature structure, allow- 
ing the temperature to drop to <1600K within a few 
photospheric radii. As it is optically thick, the CO- 
fundamental band forms at several times the photo- 
spheric radius, while the first overtone forms close to the 
photosphere (defined to be where the Thomson scatter- 
ing optical depth =1). Thus, the ratio of the fundamen- 
tal to the first overtone provides a powerful means for 
constraining both temperature and density. We use our 
near-simultaneous fundamental and first-overtone obser- 
vations to investigate conditions in the SN 2004dj ejecta. 
Only a summary of first results is given here, as a more 
detailed analysis will be presented elsewher e. 

We selected the explosion model from iChieffi et alJ 
(2003) that best matched the plateau duration of 
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Fig. 3. — Comparison of the near-IR spectra of SN 2004dj 
with those of other, roughly coeval core-collapse SNe. All spectra 
have been corrected for redshift and reddeni ng. T he SN 2002hh 
spectrum has been taken from IPozzo et alJ 120051) and has been 
smoothed for display pur poses. The SN 1987A spectrum has been 
taken from IMeikle et all <1989T) . We used E(B - V) of 0.18 for 
SN 2004dj, 2.26 for SN 2002hh, and 0.19 for SN 1987A. The flux 
scale is appropriate for SN 2004dj only. 

SN 2004dj and adjusted the 56 Ni mass to match 
the radioactive tail. The parameters of this model 
were: M — 15M Q , explosion energy = 10 51 erg, and 
M( 56 Ni)=0.03M Q . Detailed radiative transport calcula- 
tions were performed by solving the time-dependent rate 
equation s for the formati on of CO and SiO molecules, 
following iLiu et all l)1992ft . We recalculated the temper- 
ature structure for the optically thin layers under the as- 
sumption of radiative equilibrium and evolved the model 
up to day 130, when the photosphere had receded well 
into the He-enriched layers, the structure of which de- 
pends on the details of several processes e.g. mixing. 
We investigated the sensitivity to various parameters of 
CO formation, by computing an exploratory grid of 36 
models with the photospheric temperature, photospheric 
expansion velocity, and density profiles above the re- 
combination zone in the range 3000 ^ T c g(K) ^ 5500, 
2000 < v eX p(kms- 1 ) < 4000, -3 < n < -9. We find a 
best fit with T eff = 5500 K and n = 7 (Fig. |2J|, values 
which corroborate a 15M Q progenitor. 

The slope of the spectral energy distribution in the 
near-IR requires an underlying continuum to adequately 
fit the data (see Fig. If this is due to dust from 

the SN, then given the relative youth of SN 2004dj, it is 
likely to be due to pre-existing circumstellar dust, rather 
than newly condensing dust. A more likely source of this 
continu um is the underlying cluster . Using Starburst99 
models, iMafz-Apellaniz et alJ 1)20041) infer that the clus- 
ter contains «12 red supergiants, which is broadly con- 
sistent with the near-IR flux from our data, as well as 
the pre-explosion J and H magnitudes. 

3.2. Stable Nickel 

The other conspicuous feature in the SN 2004dj mid- 
IR spectra is a prominent emission feature at 6.67/im 
which we identify as the [Nill] A6.634/im fine-structure 
line. This is produced by the a 2 D 3 / 2 -a 2 -D5/2 transition 
to the ground state. That this line is already visible 
at the start of the radioactive tail, constrains the de- 
gree of mixing of the ejecta. The line is unresolved, 
indicating an expansion velocity of < 3500km s -1 . Its 
intensity at 106d is 1.1 ± 0.1 x 10 -14 ergs _1 cm~ 2 , ris- 
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Fig. 4. — Comparison between observed and model profiles for 
the vibrational bands of CO on day 130. The models are identified 
by the (T ef j ,n) tuples (see the colour version of this figure). To fit 
the observations, we have added a 1200 K blackbody component. 
Both n and T e fj have been varied to demonstrate the sensitivity 
of the emitted profiles to the physical conditions. Best agreement 
is obtained with T e g = 5500 K and n = 7 (red curve) which are 
close to those of the explosion model (see text). 



ing to 1.5 ± 0.15 x 10~ 14 ergs~ 1 cm~ 2 on day 129. We 
used this line to estimate the mass of Ni + . The criti- 
cal d ensity for the 6.63 um transition is 1.3 x 10 7 cm -3 
(e.g. iWooden et al~l 119931) . At an epoch of -100 d it 
is likely that the electron densi t y was above this value. 
For example. IClocchiatti et al~l l|1996|) estimate an elec- 
tron density of 2.5 x 10 9 cm -3 within 2500kms~ 1 for the 
Type II SN 1992H at 100 d. On this basis, we deem a 
simp le LTE treatment to be valid. A-values were taken 
from lQuinet fc L eTjoj | u3iej | r£|j[199 6|) an d partition func- 
tion values from Halcn ka et al. I l)2001f l. The Ni + mass 
was estimated for 3000 K and 6000 K. As expected, the 
result was rather insensitive to temperature. At 106 d 
and T=3000K we obtain M(Ni+)=1.7x 10~ 4 M Q . Sim- 
ilar values were obtained with T=6000K. At 129 d we 
obtain M(Ni + )=2.2x 10 -4 M©, again with similar values 
for both temperatures. Sobolev optical depths at the two 
epochs were ~0.2 and 0.4, respectively. 

Given the clear presence of the [Mil] A6.634 /xm line, 
we examined the spectra for the [Ni I] A7.50 /im line 
which is produced by the a 3 F 3 -a 3 F 4 transition to ground. 
While there is little sign of this line at 106 d, by 129 d 
there is a feature whose red wing corresponds exactly to 
the expected location of the 7.50 /xm line. We therefore 
suggest that the 7.50 /im line is present, but blended with 
another line to the blue. The most likely candidate is Pfa 



(7.46/xm). There may also be contributions to the ob- 
served feature from Hu/3 (7.50/xm) and H 7 _n (7.51/xm). 
Consequently direct measurement of the [Nil] A7.50/im 
line intensity is impractical. In a future paper, estimates 
will be made of the flux contribution from the H I lines 
in order to determine the intensity of the [Nil] line. 

At 129 d, virtually all the Ni must be made up of stable 
isotopes, dominated by 58 Ni. The derived mass can be 
regarded as a firm lower limit for the total stable Ni mass. 
As already indicated, it is likely that flux from Ni° is 
also present. Moreover, the [Nill] line is close to being 
optically thick. The presence of a significant continuum 
suggests that yet more Ni lies below the photosphere. 

3.3. Constraints on the progenitor mass 

Following lHamuv 1 l|2003j) . we use the V-band luminos- 
ity of the exponential tail of the light curve 2 after 100 d 
to estimate the mass of 56 Ni. We find a mean mass of 
~0.022 My which suggests a progenitor mass of > 10 M© 
l)Wooslev fc Timmes Ill99"fi) and is consistent with the 

v alue used in A3. II 
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1)19961) predict stable Ni masses for 
core-collapse SNe having progenitor masses of 13-25 M Q . 
Up to 20 M Q , they predict masses of 0.007-0.013 M Q , but 
a much lower value of 0.002 M© for a 25 M Q star. At this 
early epoch, our lower limit does not yet constrain this 
mass range. However, as the SN expands, we will be able 
to to obtain a more definitive Ni mass estimate. 

By fit ting a variety of cluster spe ctral energy distri- 
butions, iMaiz-A pcllani z'et alJ l)2004D have inferred that 
Sandage 96 is a 13.6 Myr cluster with a turn-off mass of 
15Mq. Although they favour a red supergiant progenitor 
for SN 2004dj , they could not completely r ule o ut a blue 
supergiant progenitor. Our CO analysis ( fl3.1|l supports 
a red supergiant progenitor. This is indirectly supported 
by the lack of PAH features (at 6.2, 7.7, 8.6, and 11.3 /im) 
which trace the far-UV stell ar flux and therefor e, the 
young, hot stellar population l)Peeters et alJ l2004). 

This work is based on observations made with the Spitzer 
Space Telescope, which is operated by the Jet Propul- 
sion Laboratory, California Institute of Technology under 
NASA contract 1407. Support for this work was provided 
by NASA through an award issued by JPL/Caltech. 
R.K. acknowledges support from the EC Programme 
'The Physics of Type la SNe' (HPRN-CT-2002-00303). 
We thank C. Gerardy for useful discussions, and M. 
Pozzo for providing the SN 2002hh spectrum (Fig. |3J). 
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